Report for Model 6

— for both jointly and separately

2024-10-08

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.03       1.03
## dh0           1.18       1.64
## dh1           2.99       6.93
## dh2           4.48       9.72
## dl0           1.11       1.43
## dl1           2.17       4.55
## dl2           3.02       6.23
## dl3           1.15       1.55
## 
## Multivariate psrf
## 
## 2.93

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1227.54301 23434.02
DIC3 1162.80834 20673.75
PWAIC 46.27266 1002.20
WAIC 1192.59910 21062.58

2. Separate Modeling of High-Level

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.00       1.02
## dh0           1.23       1.79
## dh1           1.12       1.40
## dh2           1.02       1.09
## dl0           1.11       1.12
## dl1           1.15       1.51
## dl2           1.03       1.11
## dl3           1.19       1.37
## 
## Multivariate psrf
## 
## 1.21

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

H2
DIC 1377.52079
DIC3 1247.01158
PWAIC 90.84391
WAIC 1316.35844

3. Separate Modeling for Low-level

3.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

3.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## Devsum        1.01       1.03
## LDevsum       1.00       1.02
## dl0           1.11       1.12
## dl1           1.15       1.51
## dl2           1.03       1.11
## dl3           1.19       1.37
## 
## Multivariate psrf
## 
## 1.08

3.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

3.4 WAIC results

L7
DIC 23110.679
DIC3 20729.655
PWAIC 1040.447
WAIC 21133.344